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^\ ■ We study the exchange of kinetic energy between translational and rotational degrees of freedom 

' for inelastic collisions of rough spheres. Even if equipartition holds in the initial state it is immedi- 

, ately destroyed by collisions. The simplest generalisation of the homogeneous cooling state allows 

^ ■ for two temperatures, characterizing translational and rotational degrees of freedom separately. For 

times larger than a crossover frequency, which is determined by the Enskog frequency and the initial 
temperature, both energies decay algebraically like with a fixed ratio of amplitudes, different 
from one. 
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,—1 ' Kinetic theory of inelastically colliding particles has become a subject of growing research activity, motivated partly 
,.Ch ' by renewed interest in granular materials. A Boltzmann equation and the Enskog variant of it have been formulated 
for inelastically colliding particles with normal and tangential restitution . Whereas the derivation of the kinetic 
equation is relatively straightforward, the methods of solution which were developped for elastic collisions cannot be 
taken over to the inelastic case because there is no simple stationary or local equilibrium distribution, around which 
one could expand the nonlinear kinetic equation. One simple distribution is the homogeneous cooling state (HCS) 
which depends on time only implicitly via the average kinetic energy T{t). The latter is predicted to decay like 
c/2 ^ for large times. The homogeneity assumption is certainly violated when clustering occurs and, in particular, if 

inelastic collapse happens. Nevertheless we have found recently ||] in a model of inelastically colliding rods that the 
kinetic energy follows on average a behaviour, even if clustering occurs. This suggests that the above scaling law 
may be a useful approximate description, even when the assumptions of HCS break down. 
' In this note we show that rotational and translational energy relax with different rates in general. Once friction is 
(— I included, HCS with a single time dependent temperature is no longer consistent with the time evolution of translational 
Q and rotational energy separately. Instead one has to introduce two temperatures which characterize translational and 
O \ rotational degress of freedom separately. Both are found to fall off like with the ratio approaching a constant 

value, which is determined by the coefficients of normal and tangential restitution, 
^f-^ \ We briefly recall the collision dynamics of hard spheres with normal and tangential restitution. These results can 
^ . for example be found in Cerginani Q. We consider two spheres of equal diameter a, mass M and moment of inertia 
CO ' /. The unit-vector from the center of the first sphere to the center of the second is denoted by n and velocities and 
angular velocities before collision by Vi, t>2, tc^i and uj2- The relative velocity of the contactpoint before collision is 
^5 ' given by V = i;2 + %n x — I'l + %n x oji . Normal and tangential restitution determine the relative velocity after 
, collision according to 

; n • y = -e(n ■ y) with ee[0, 1], (1) 

n X y' = -/3(n X V) with {3 G [-1, 1]. (2) 



Using the property of conserved linear and angular momenta one obtains for the velocities after collision 



■ v^=vi- i]tvi2 - (?7„ - r7t)(n • Vi2)n + rjt^h x (wi +^2), 

O' ' , ^ . ^ a ^ . . 

O . ^2 = f 2 + VtVi2 + (Vn - Vt)in ■ vi2)n ~ r]t-n x (wi +^2), 



iv^ = LJi -ritn X V12 + —n x (n x (lji + u;2)), 

OjK fC 

U}2= ^2 T^*" ^ ''^12 + "T" X ('^^i + '^2))- (3) 



X 

^, U.2 - U.2 - — ^ ^12 -r - 

with V12 — Vi — V2 and parameters k := (k=0.4 for homogeneous spheres), ry„ := i^^, and rjt := 2[i+k] • 

We consider a system of N classical particles, confined to a 3-dimensional volume V and interacting via a hard 
core potential. Each particle is characterized by its position ri(t), its linear momentum Pf{t) = Mvi{t) and angular 
velocity uJi{t). The time development of a dynamical variable A = A{{ri{t),p^{t),uJi{t)}) is determined by a Pseudo 
Liouville-operator 

A{{n,p^,LO,},t) = exp{tC+t)A{{r,,p„uj,},0) for t > 0. (4) 
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Such a Pseudo Liouville-operator was first introduced for hard (perfectly smooth) spheres by Ernst et al. |10| and 
subsequently applied to the calculation of transport coefficients of a hard sphere fluid (see e.g. [0 and references 
therein). Noije and Ernst have generalised the formalism to inelastic collisions with normal restitution. Here we 
extend these results to rough spheres with normal and tangential restitution. 

The Pseudo Liouville-operators £+ consists of two parts — Co + Cj^. The first one, £o, describes free streaming 
of particles 

^o = -]^EP"-^'-" (5) 

n 

and the second one, = ^ ^n=^m T+{nra), describes hard-core collisions of two particles by 



T+{nm) = -^(p„,„ • f-„„)e(-p„„ • f„™)(5(|r„„| - a)(6+„ - 1). (6) 

The operator 6+„j replaces the linear and angular momenta of two particles n and m before collision by the correspond- 
ing ones after collision. Q{x) is the Heaviside step-function and we have introduced the notation p„„j = p„ — p„j, 
rnm = rn- r,n and r = r/\r\. 

Eq. (^) has a simple interpretation. The factor (p„„j • f„,„)/Af gives the flux of incoming particles. The 0- and 
(5-functions specify the conditions for a collision to take place. A collision between particles n and m happens only 
if the two particles are approaching each other, which is taken into account by 0(— p„„ • rnm)- At the instant of a 
collision the distance between the two particles has to vanish, as expressed by i5(|r„TO| — a). Finally (6„m~l) generates 
the change of linear and angular momenta . 

The ensemble average of a dynamical variable is defined by 

< A >t= I dTp[Q)A[t) = fdTp{t)A{0) - f Y[{dr,dp^duj,)Y[e{\nj\ - a)pit)A{0). (7) 

Here p{t) ~ exp {—iC\^_t) p{0) is the iV-particle distribution function, whose time evolution is governed by the adjoint 

C\_ of the time evolution operator The quantities of interest are the translational and rotational energies per 
particle 



E, 



ot 



as well as the total kinetic energy E = En + i^rot ■ 
As a first step we compute the initial decay rates 

^ < Etr >t=o = < iC+Etr >t=o= (10) 
d 

— < Ei-ot >t=0 — < i^+Erot >t=0— i'rot (H) 

assuming that the system has been prepared in a thermal equilibrium state p(0) oc exp(— i?/T) which is characterized 
by the temperature 3T/2 ^< E >t=o. The change of translational and rotational energy are given by 



(7rM)i/2 \l + k 

Here g(a) denotes the pair correlation at contact and uq = N/V. Rotational energy is conserved in two cases, for 
either perfectly smooth spheres (/3 = —1) or perfectly rough spheres (/3 = -1-1). Translational energy is only conserved 
if in addition e = 1. For all other values of the parameters e and (3 the translational and rotational energy decrease 
linearly with time but with different rates. This implies that after a small time interval At equipartition among 
rotational and translational degrees of freedom no longer holds 
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< Etr >At = < Et, >Q -VtrM 
< Erot >At ~ < Erot >0 — ^'rotAt. 



(14) 



Hence for collisions with friction {(3 ^ —1) the homogeneous coohng state is not consistent with the time evolution of 
translational and rotational energy separately. To generalise the concept of a homogeneous cooling state to collisions 
with friction we introduce two temperatures rtr(i) = 2/3 < i?tr >t and T'rot(^) — 2/3 < i?rot >t- We keep the 
assumption of spatial homogeneity and assume that both linear and angular momenta are normally distributed with 
in general different time dependent widths or temperatures 



1 ^ 

p{t) cx exp -- ^(Mi;2/7^tr(i) + 



(15) 



The above distribution allows for a calculation of the dacay rates of Ttr(i) and Ti.ot(i) for arbitrary times, resulting 
in two coupled differential equations: 



3/2 



dt 

d rr. 7 2^3/2 Vt 



dt 



(-1) 



Jrot-ttr 



7 2rjnl/2rjn 



(16) 
(17) 



with 7 = {16/3)a'^nog{a)^/TrjM. 

As we have seen above the energies decrease linearly for short times. For large times both, translational and 
rotational energy fall off algebraically like with however different amplitudes. The amplitudes can be determined 
analytically, e.g. by solving the differential equation for 



dTrr. 



-Ttr ((1 - e")/4 + rjtjl - Vt)) + TrotyVk 
TtrVt/k - Troti^ - Vt/k)vt/k 



(18) 



This equation is solved by a constant ratio Ttr/Trot = —Ci/c2 + \/l + (ci/c2)^ with ci, C2 given in terms of e, f3 and k 
by ci = (1 -e2)/4+ (1 -/32)(fc- l)/(4fc + 4) and C2 = 2fc(l + /3)2/(2 + 2fc)2 in agreement with ref. [||. The crossover 
time to between the linear regime and the algebraic decay is determined b y 7 and Ttr(O) according to ^ oc 7\/Ttr(0). 
The full solution has been obtained by numerical integration of eqs.( p6|Jl7| ) and is shown in Fig. (|l|) for k = 0.4, 
7 = 1, e = 0.4, (3 = —0.6 and Ttr(O) = ri-ot(O) = 20. The asymptotics is indicated by a dashed-dotted line and the 
inset is a blow up of the short time behaviour. 




FIG. 1. Time decay of the translational and rotational energy. Parameters were chosen according to: k = 0.4, e = 0.4, 
P — —0.6 and rtr(O) — rrot(O) = 20. We have introduced a dimensionless time argument r — tjy/Ttv{0). The dashed-dotted 
line indicates the asymptotics 
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The deviation of the ratio of rotational and translational energy from 1, i.e. 1 — Tiot/Ttr, is shown in Fig. 
two different values of (3 and e. 
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FIG. 2. Ratio 1 - Trot/Ttr as a function of time with k = 0.4, 7 = 1, Ttr(O) = rrot(O) = 100, e = 0.6, (3 = -0.4 and e = 0.4, 
/3 — —0.6. Unit of time was chosen as in Fig. (|l]). 

Cooling of a granular gas has been investigated by several groups ||l6|,0 modelling collisions with a normal coefficient 
of restitution. The behaviour of the total energy has been confirmed for a wide range of parameters such that the 
system is stable to shear fluctuations |l7) and inelastic collapse. Luding ||l8| has considered a more detailed model of 
inelastic collisions including tangential restitution and Coulomb's law of friction. It is straightforward to generalise 
the model to include Coulomb friction, one just has to modify the updating rules (eq. (|^)) accordingly. Driving the 
system with a vibrating wall is an open problem so far. One may also generalise the approach to nonspherical objects, 
e.g. hard rods or needles ||l9|, for which elastic collisions have been analysed by Frenkel and Maguire pc| ]. 

We wish to thank Timo Aspelmeier and Kurt Broderix for useful discussions and one of us (M.H.) the Land 
Niedersachsen for granting a scholarship. 

Note added in proof: After this paper was completed we learnt about related work by S. McNamara and S. Luding 
who also study the ratio of translational and rotational energy in the HCS. 
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